Maxwell Model of Traffic Flows 
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We investigate traffic flows using the kinetic Boltzmann equations witfi a Maxwell collision 
integral. This approach allows analytical determination of the transient behavior and the size 
distributions. The relaxation of the car and cluster velocity distributions towards steady state is 
characterized by a wide range of velocity dependent relaxation scales, R^^^ < r(«) < R, with R 
the ratio of the passing and the collision rates. Furthermore, these relaxation time scales decrease 
with the velocity, with the smallest scale corresponding to the decay of the overall density. The 
steady state cluster size distribution follows an unusual scaling form Pm ~ (m)~^'I'(m/{m)^). This 



distribution is primarily algebraic, Pm ~ m for m <^ 
PACS numbers; 02.50-r, 05.40. +j, 89.40+k, 05.20.Dd 



(m) , and is exponential otherwise. 



I. INTRODUCTION 



Traffic flows exhibit a variety of collective behaviors 
typical to non-equilibrium systems The observed 

phenomenology is rich and includes shock waves, traffic 
jams, clustering, and synchronized flow [^-Q. A num- 
ber of models and theoretical approaches including fluid 
mechanics] |,^D , cellular automata , particle hop- 

pingjl8-20|, ballistic motion ||2^-|2^, and optimal veloc- 
ity p9|-|31| are used to describe these observations. Yet, 
different approaches have different virtues, e.g., kinetic 
theory is more appropriate for dilute flows, while fluid 
mechanics is more appropriate for dense flows. 

Here, we focus on the kinetic description of traffic. 
Previously, we introduced a microscopic ballistic motion 
model and used it to derive Boltzmann Equations (BE) 
for traffic on a one lane roadway ||2^ . A generalization to 
situations where passing is allowed shows that a transi- 
tion from a low-density "laminar" flow to a high-density 
"congested" flow generally occurs p5| , p6t . The resulting 
BE are technically difficult, and a number of important 
questions remain unresolved including the transient char- 
acteristics and the cluster-size distribution. Indeed, pre- 
vious studies addressed only steady state properties and 
the results concerned mainly the velocity distributions. 

Our goal is to obtain these relevant flow characteris- 
tics. To this end, we propose an approach inspired by 
Maxwell's classical model, widely used in kinetic theory 
|^,|3|. This Maxwell approach uses a velocity indepen- 
dent collision rate, thereby considerably simplifying the 
analysis. In fact, upon transforming the kinetic equa- 
tions from integral into differential ones, the Maxwell 
model results in first order differential equations while 
the Boltzmann approach leads to second order equations. 

We will show that the Maxwell approximation is faith- 
ful to the nature of the original traffic equations as it 
qualitatively reproduces transient characteristics for no 
passing zones, as well as steady state characteristics for 
passing zones. We further find that the cluster veloc- 
ity distribution approaches its steady state according to 
a wide spectrum of relaxation scales, with the smallest 



describing decay of the overall cluster density. Further- 
more, the size distribution is characterized by a strong 
algebraic tail for small and average sizes, while it is ex- 
ponentially small for large size. 



II. THE MAXWELL MODEL 

The ballistic motion approach models the basic pro- 
cesses underlying one lane traffic ffows: passing and slow- 
ing down due to clustering. The main assumption is that 
each driver has a fixed intrinsic velocity. The driving 
rules are as follows: A car moves with constant intrinsic 
velocity on a one lane road until it overtakes a smaller 
velocity car or a cluster. After such an encounter, or "col- 
lision" , the incident car immediately adopts this smaller 
velocity, thereby joining a cluster. Cars may also resume 
driving with their intrinsic velocities, and such passing 
events are assumed to occur with a constant rate. This 
model is an idealized description of one lane traffic ffows 
as several time and length scales including the actual 
collision time, the passing time, and the car size are ne- 
glected. 

Let P{v,t) be the density of clusters moving with ve- 
locity V at time t, and let Po{v) be the intrinsic veloc- 
ity distribution. Natural initial conditions where cars are 
randomly distributed in space and drive with their intrin- 
sic velocities are imposed, i.e., P{v,0) — Po{v). Under 
the assumption that space and velocity remain uncorre- 
lated, a mean-field Boltzmann equation is written 



dPiv,t) . 



dt 



t^'[P,,{v)-P{v,t)] 

-P{v,t) / dv'U{v,v')P{v' ,t). 
Jq 



(1) 



The first term on the right-hand side represents cars es- 
caping their respective clusters with a constant rate t^^. 
The last term accounts for decrease in the cluster den- 
sity due to collisions. For traffic ffows the collision rate 
should read U{v,v') = v — v' . For such a collision rate, 
steady state properties of the velocity distributions can 
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be obtained by transforming the rate equation into a sec- 
ond order nonlinear differential equation |25| . However, 
a number of important characteristics including the size 
distribution and time dependent properties appear to be 
analytically intractable. 

We propose using a constant collision rate, U{v, v') = 
Mq to simplify the above traffic equations. Similar 
in spirit approximations, termed the Maxwell Model 
(MM), proved useful in kinetic theory [^,|3j. A nat- 
ural choice for the constant rate uq is the average ve- 
locity difference, uq ~ ^ 'v') (v). One may won- 
der whether such an approximation is reasonable for 
traffic flows. Ignoring P{v',t) in the collision integral 
I{v) = /p dv'U{v,v')P{v' ,t), we have I{v) cx v for the 
MM, while I{v) oc for the BE. Hence, the integral 
remains an increasing function of the velocity. Further- 
more, cars must slow down before a collision, and there- 
fore, the collision rate should be slower than linear. The 
MM can actually be considered as the limiting case of 
zero deceleration, while the linear rate corresponds to 
the limit of infinite deceleration. 

Let Co be the initial car concentration, vq the aver- 
age intrinsic velocity, Iq^ the passing rate, and uq the 
collision rate. Introducing the dimensionless velocity 
v/vq —> V, space xcq and time coVot — > t variables 

normalizes the initial concentration and typical velocity 
to unity. The master equation (|I|) is characterized by two 
dimensionless numbers 

The normalized collision rate, v — uq/vq, merely rescales 
time. Thus, it is set to unity without loss of generality. 
The number R = cqWo^o = ^osc/^coh equals the ratio of 
the two elementary time scales, the escape time tcsc = to, 
and the collision time icoi = (cqKo)"^- This number, 
termed the "collision number" , plays an important role 
in determining the nature of traffic fiows. 

We will show below that the Maxwell model is com- 
pletely solvable. Although quantitative results of the MM 
may differ from the BE, they faithfully reproduce the 
qualitative behavior of the traffic equations. 



III. THE CLUSTER VELOCITY DISTRIBUTION 

We start with steady-state and time dependent proper- 
ties of the cluster velocity distribution. Let us introduce 
the auxiliary function 



Q{v,t) = R- 



dv'P{v',t), 



(3) 



which gives the cluster distribution via differentiation 
P{v,t) = dQ{v,t)/dv. This auxiliary function enables 
us to reduce the integro-differential Eq. (j^) into the dif- 
ferential equation 



d_dQ ^ 1 dQo 
dt dv R dv 



ov 



(4) 



This equation can be integrated over v, and using the 
boundary condition Q\v=q = R^^ we find 



dQiv,t) Qoiv) Q^v.t) 



dt 



R 



1 

2R^' 



(5) 



Integrating Eq. the auxiliary function is obtained 
explicitly for arbitrary initial conditions 



Qiv,t) 



1 + A(w)e-*'3<»W 
1- A(w)e-*'5~W 



(6) 



with notation = [Qoiv) - Q^{v)]/[Qq{v) + Qodv)]. 
Here we use the subscript oo to denote steady state. The 
steady state auxiliary function Qoo{v) = Q{v,t = oo) is 
given by 



Qoo{v) = R 



2R I dv'Poiv') 

'0 



1/2 



(7) 



Since the concentration is obtained from Q{v,t) using 
c(t) — \im..^,-^oo[Qiv,t) — R^^], and since the cluster ve- 
locity distribution is obtained by differentiation, Eq. (|^) 
represents a complete explicit solution of the Maxwell 
model. 

We first examine steady state properties of the cluster 
velocity distribution. Comparing with the correspond- 
ing behavior emerging from the BE will allow us to test 
the utility of the Maxwell model. Evaluating the infinite 
velocity limit of the auxiliary function gives the overall 
cluster density 



R 



2R- 1 



(8) 



A remarkable feature of the steady state cluster density 
is that it is a function of the collision number only. Such 
independence of the initial velocity distribution has been 
observed in a few other ballistic aggregation problems 
|9|,|§. Note that = l-R/2 + 0{R^) for i? < 1, i.e., 
the difference from the initial density is of order R in the 
laminar flow regime. In this study, we will focus on the 
complementary nontrivial limit of congested flows, i.e., 
R ^ 1. Here, the cluster concentration is significantly 
reduced, Coo ~ R~^^^, and large clusters with an average 
size (m) = c -1 ~ form in agreement with the BE 

results. 

The cluster velocity distribution is obtained from 
Eq. M) by differentiation 



Pooiv) = Poiv) 



1 + 2R dv'Poiv') 



-1/2 



(9) 



When R 1, the difference between the initial and the 
steady state distributions is of order R. This indicates a 
laminar flow regime when the correction due to collisions 
is small and can be obtained by expanding the solution 
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perturbatively around the initial state. When i? ^ 1, 
we use the notation Io{v) = dv' Po{v') and write the 
leading behavior of Eq. (0) as 



Po{v) 



(10) 



The two limiting behaviors match at the threshold veloc- 
ity V* which is found from 1 ~ RIo{v*) = R dv Po{v). 
In agreement with the Boltzmann approach [ p5p^ ], a 
boundary layer structure is found for the velocity distri- 
bution, where in the inner region the original distribu- 
tion prevails, while in the outer region, the distribution 
is substantially reduced. The average cluster velocity re- 
mains of order unity. Additionally, the suppression of 
the fastest velocities is proportional to the concentration, 
again in agreement with the BE results. We conclude 
that although the MM differs quantitatively from the ex- 
act BE behavior, it qualitatively reproduces the steady 
state behavior. 

We turn now to analyzing the transient properties and 
in particular the approach towards steady state. The 
time dependent concentration reads 



C(t) = Co 



1 + 56-*/^"= 
' 1 - Se-*/^c ' 



(11) 



with the constant B = A{oo) = (1 — Coo)/(l + Coo) and 
the relaxation time Tc = R/ \/2R + 1 corresponding to 
the concentration decay. We see that the cluster concen- 
tration exponentially approaches its steady state value 



c{t) 



1 + 2Be~*/^'^ 



(12) 



As the distribution changes slightly in the laminar phase, 
relaxation times remain of order unity when R 1. 
However, for congested flows the relaxation time diverges 
with the collision number Tc ~ R^^^. 

The explicit time dependent auxiliary function allows 
determination of relaxation properties of the cluster ve- 
locity distribution. In the long time limit Eq. (|^) reads 



Qiv,t)^Qooiv) [l + 2A(t;)e-*/^('') 



(13) 



with the velocity dependent relaxation time t{v) — 
l/Qoo(w)- Thus, the steady state properties are re- 
flected in the transient characteristics. The velocity de- 
pendence of the relaxation time t(v) becomes especially 
pronounced for large collision numbers where it exhibits 
the following boundary layer structure 



f . ( R w < w*; 

^^""'^ \ [R/Io{v)]'/^ v^v*. 



(14) 



For sufficiently small velocities, the collision integral 
is negligible, and the relaxation time R, suggested by 
Eq. (g) holds. While small velocities are governed by 
(almost) velocity independent relaxation scales, large ve- 
locities are characterized by velocity dependent decay 



rates. Furthermore, a large range of relaxation scales 
exists i?^/^ < T < R with the larger relaxation scales 
corresponding to smaller velocities. This is consistent 
with dimensional arguments that time and velocity are 
inversely related. Interestingly, the smallest possible re- 
laxation scale corresponds to the overall cluster density. 

One anticipates that the relaxation time r(u) also gov- 
erns the decay of the cluster velocity distribution P{v,t). 
This is indeed the case. To obtain explicit expressions 
we first simplify Eq. (O), 



Qiv,t) 



(15) 



Differentiating with respect to v gives the cluster velocity 



2RPo{v)Io{v)e-*/^^''^ 
P{v,t)-P^ ~ { -RPo{v)liiv)[te-'/^(-^] 
-2R-^Po{v)[te''/'-^'''>] 



V < V**; 

V > V*; 



with t{v) given by Eq. (14). The expressions match 
at the boundary velocities which are determined from 
RIo{v*) ^ 1 and tIo{v**) ~ 1. Only for velocities slower 
the decaying boundary velocity v**(t) is the correction to 
the cluster density positive. This is surprising since both 
the overall cluster density and the auxiliary function ex- 
hibit positive corrections, as one would naively expect 
since Po{v) > Poo{v). 

Since the relaxation times diverge with increasing i?, 
an intermediate behavior should emerge in the time range 
t <C V^. In this regime, the system has not yet "real- 
ized" that passing is allowed, and the behavior should 
agree with the no passing case where i? = oo. By di- 
rectly solving Eq. (|^) with R~^ = one finds 



Piv,t) 



Po{v) 



[i + ^tio{v)y 



(16) 



For arbitrary intrinsic velocity distribution, a scaling 
asymptotic behavior emerges 



(vit)) \{v{t))J' 



(17) 



with the cluster concentration c{t) ^ t^^ and the average 
velocity determined by Hq {{v{t))) ^ 1. We see that the 
average velocity in the no passing case is proportional 
to the time dependent boundary velocity: {v{t)) ~ v**. 
When the leading small velocity behavior of the intrinsic 
velocity distribution is algebraic, Pq{v) ~ as u ^ 0, 
the average velocity decays as a power law in time, 
{v{t)) ~ with f3 = l/{fi + 1). Comparing with the 
exact behavior in the no passing limit of ballistic motion 
model, we see that the overall scaling picture is repro- 
duced, while the quantitative details and in particular 
the scaling exponents are different. 

To summarize, explicit expressions for all cluster prop- 
erties are possible in the realm of the MM. The relaxation 
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towards steady state occurs in two stages. The early 
one, t <C VR, corresponds to a no passing intermediate 
asymptotics. Then, the passing mechanism comes into 
play, and the system approaches steady state. This lat- 
ter relaxation is nontrivial in several ways. The decay is 
non-uniform in time as a wide range of time scales are 
observed. It is also non-uniform in velocity as the cluster 
velocity distribution involves three layers of greatly dif- 
ferent width, i.e. it exhibits the so-called "triple deck" 
structure [Q. The first layer v <C v**{t) (referred to 
as the lower deck) shrinks with time and the velocity 
distribution in this deck approaches the steady state ex- 
ponentially from above. In the middle and upper decks, 
the approach towards steady state is from below and has 
a linear in time correction to the exponential decay. 



IV. THE CAR VELOCITY DISTRIBUTION 

The cluster velocity distribution does not provide the 
observed distribution of car velocities since all clusters - 
large and small - are taken with equal weight. In what 
follows, we concentrate on the car velocity distribution, 
which determines basic properties such as the flux. 

Within the framework of the MM, the car velocity dis- 
tribution satisfies 



dG{v,t) 

di 



R-^[Po{v) ^ G{v,t)] 



(18) 



-G{v,t) dv'P{v\t)+P{v,t) dv'G{v',t). 

Jo Jv 

The escape term is the sum of a gain term R^^[Po — P] 
and a loss term —R^^[G — P]. In a collision between two 
clusters, all cars belonging to the faster cluster acquire 
the slower cluster velocity. Thus, in both collision terms 
the argument of P is smaller than the argument of G. In 
contrast with Eq. (H) collisions can now lead to a gain in 
the car velocity distribution. Although the integration 
limits resemble those of the previous kinetic equations 
p3[ , the collision terms are different, a reflection of the 
different treatment of cars and clusters in this theory. 
One can verify that Eq. (|l^) conserves the car density 
1 ~ dv G{v,t). An alternative approach for obtain- 
ing G{v, t) involves the conditional velocity distribution 
P{v,v',t). This more detailed distribution can also be 
used to verify G(v,t), and for completeness, we detail its 
derivation in Appendix A. 

Let us introduce the auxiliary function 



dv'G{v',t). 



(19) 



In terms of the auxiliary functions g, Q, and Qo, Eq- ( |l^ 
becomes 



d dg{v,t) d 
dt dv dv 



R 



(20) 



Integrating over the velocity and using gQ{v) = Qo{oo) - 
Qo{v) gives the master equation 



dgiv,t) 
dt 



R-^go{v)~g{v,t)Q{v,t). 



(21) 



We first analyze the steady state properties which are 
obtained immediately from Eq. (|2l|) 



9o{v) 

RQoo{v)' 



(22) 



Interestingly, this auxiliary function and the cluster ve- 
locity distribution experience the same relative reduc- 
tion at the steady state, goo{v) / go{v) = Poc{v)/Po{v) = 
1/RQqc{v)- Differentiating goo{v), we get 



G^iv)^Po{v) 



l + R + RIo{v) 
[l + 2RMv)] 



3/2 ■ 



(23) 



In the congested phase, i? ^ 1, the car velocity distribu- 
tion has the following limiting behaviors: 



Goo(«) 



RPoiv), 



V < V*: 



R-^/^Po{v)I-^/\v), v^v* 



(24) 



Thus, while the fast tail decay agrees with the 

Boltzmann equation approach [ p5[ , the slow tail enhance- 
ment R is larger than the Boltzmann result where this 
enhancement is of the order R" with < a < 1. 

The car velocity distribution immediately gives the av- 
erage size of a i;-cluster 



l + R + RIojv) 
l + 2RIo{v) 



(25) 



obtained from {m{v)) — G{v)/P{v). As expected, the 
average cluster size decreases with the velocity. The av- 
erage cluster size obeys the bounds 1 < {m{v)) < I + R, 
with the upper (lower) bound achieved by the slowest 
(fastest) clusters. An additional quantity, immediately 
derived from the car velocity distribution is the ffux, 
Joo — J dvv Goo (v) ■ One can use Eq. ( |2^ ) to find 



dv ■ 



1 - Io{v) 



+ 2RIo{v) 



(26) 



In the congested phase, the ffux is proportional to the 
threshold velocity, Joo ^ v* , in agreement with the Boltz- 
mann equation results. 

We now focus on the time dependent behavior. Inte- 
gration of equation ( ^l|) gives an explicit expression for 
g{v, t) (for a derivation, see Appendix B) 



g{v,t) Q{v,t) 



9oo{v) 



Qoo{v) 



1 



loiv) 



(27) 



The relaxation of g follows directly from the relaxation of 
Q since g{v, t) — goo{v) o: Q{v, t) — Qoo(y) when t oo. 
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Using Eq. (15), we evaluate the leading relaxation behav- 
ior of g{v, t) 

f 2i?go(f)/o(w)e-*/"(") v<^v**; 

[ -2i?-i5o(«)[ie-*/^('')] V > V*. 

Interestingly, the relaxation of the (properly normalized) 
cluster and auxiliary car velocity distribution are iden- 
tical,^ [g{v,t) - goo]/goiv) ^ [Piv,t)^ - Po.]/Po{v). Re- 
laxation of the car velocity distribution is obtained from 

G = —dg/dv 



G{v,t) 
GUv) 



-/2(z;)[t2e-*/-Ml 
-/o(w)i?-H^'e-*/"(")] 



V > V*. 



Thus an exponential relaxation with a velocity depen- 
dent time scale t{v) underlies the approach of all velocity 
distributions towards steady state. The car velocity dis- 
tribution exhibits the triple deck structure similar to that 
of the cluster velocity distribution. Some details change 
however; for example, in the middle and upper decks the 
prefactor further slows down the decay of G{v, t). The 
car velocity distribution approaches its steady state al- 
ways from below. 



V. THE SIZE DISTRIBUTION 

An important characteristic of traffic flows, the cluster 
size distribution, has been determined analytically only 
in the no passing limit ||2^. We now address this issue 
in the framework of the MM. Let us consider Pm{t) the 
cluster size distribution which evolves according to the 
following rate equation 



dPrnjt) 

dt 



R-\mP,^+i{v,t) - (m - l)P„(t)] (28) 



+ R-Hm.l[l-c{t)]-c{t)Pra{t) + ]^ J2 ^^W^jW- 

Terms proportional to account for escape, while 

the rest represent collisions. The overall collision rate 
experienced by a cluster, c{t), is velocity-independent. 
These rate equations are reminiscent of an aggregation- 
fragmentation process [^j37|. Indeed, collisions lead to 
cluster aggregation while passing events split clusters. 

Since aggregation and fragmentation are opposite 
mechanisms, their combined effect generally leads to a 
steady state. We leave the ambitious task of a complete 
solution for the future, and restrict our attention to the 
steady state where Eqs. (ES) read 



jP„i = R ^[mPni+i - (m - l)Pr„] 

+ i?-i(5„,,i(l-coo) + i P'P^- 

i-\-j—7n 



(29) 



It is useful to introduce the generating function 

m 

At the steady state, it satisfies the Riccati equation 



d (T 



-^z(l-z)- - =0. (31) 

1 — Coo az \ z 



j)/(i?c^) = 1 was used in obtaining 



The identity (1 - 
this equation. 

Although we could not solve these equations gener- 
ally, most of the interesting features can be obtained by 
carefully analyzing the leading terms in R. The asymp- 
totic relation Cqo — ■\/2/i? suggests that the last term 
in Eq. ( ^ ) is negligible. Solving the resulting quadratic 
equation subject to the boundary condition = 1 

gives !F = 1 — VI ~ -Z- Expanding this expression in 
powers of z we arrive at 



Prr 



2r(i)r(m + i) 



(32) 



which simplifies to Pm — (27ri?)^^/^m^'^/^ for m ^ 1. 
However, this solution does not apply for very large m, 
or equivalently near z = 1. This follows e.g. from the 
conservation of the car densit y, "mPm — 1, which 
implies that a crossover from (|32| ) to the tail behavior 
should occur at the cutoff size nic ~ (m)^ ~ R. 

To investigate the very large m behavior we have to re- 
turn to Eq. (|3l|). Fortunately, in the proximity of 2: = 1, 
i.e., when 1 — z ~ R~^, the generating function depends 
on a single scaling variable 



l-.F = Coo$(C), with C 



1 - z 



(33) 



This can be seen by balancing the leading terms in 
Eq. (31). The scaling hmction <I> satisfies the Riccati 
equation 



C<i>'(C) = C - 



(34) 



subject to the boundary condition $(0) = 0. Using the 
transformation $(C) = 4'{0 1 4'' {0 reduces Eq. ( p^ to a 
second order linear differential equation 



(35) 



This is a solvable one-dimensional Shrodinger equation 
for a particle with zero energy in a repulsive Coulomb 
potential. Indeed, a solution is found by reducing 
Eq. (|35| ) to the Bessel equation. Choosing the solu- 
tion which satisfies the appropriate boundary conditions, 
= 0,</.'(C) = 1 at C = 0, one finds 0(C) = VQh{2^) 
with Ii{x) the modified Bessel function of order one. Re- 
turning to $(C) we obtain 
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$(C) = 2C 



The last expression is derived using the identities /( {x) = 
Io{x) - x-^hix) and /^(x) = h{x) |§. 

The function <I>(C) is the Laplace transform of the prop- 
erly scaled size distribution. Indeed, Eq. (p3|) implies 



E^n^(l-^'")=C^<i>[c- 

the scaling form Pm{R) - 



■(1 — z)] whose inversion yields 



^^{c^m). Therefore, in the 
large R limit the size distribution follows the scaling form 



1 



(37) 



with (m) = 1/coo ^ \jR/'2.. The scaling function vI/(M) 
obeys $(() = /o°° dM *(M) [l - 6"^^^] . Differentiating 
both sides with respect to C, shows that <&'(C) is simply 
the Laplace transform of Af4'(Af) 



$'(C) 



(38) 



Consequently, the asymptotic behavior of the size distri- 
bution can be determined from the corresponding asymp- 
totics of $(C)- The latter are found from Eq. (p6[): 



*(C) 



/C(C 
l^/C 



-cr 



(39) 



The algebraic behavior of $(C) at large ^ implies an al- 
gebraic behavior of 5* (A/) at small M; similarly, the pole 
at C = -C* (C* = 1.445796 Q) implies an exponential 
decay for large M: 



^'(M) 



(47r)-i/2M-3/2, M<1; 



C*exp(-C*M), Af>L 
In terms of the original variables we have 



(40) 



" - \ 4C*-R-' exp{~2Cm/R), m > i?. ^ ^ 

These two limiting behaviors match ai m ^ R where 
Pm ^ R^^ ■ Additionally, the value of the cutoff size, 
•nic ~ i?, agree with our previous findings. 

In conclusion, the Maxwell equation approach allows 
explicit calculations of the size distribution. It decays al- 
gebraically with size for small and average clusters, and 
exponentially for very large clusters. The interesting as- 
pect of the size distribution concerns its scaling form. 
If the typical and the average size would be the same, 
a naive scaling argument m/{m) would under ly the size 
distribution. However, a different picture emerges where 
the scahng variable is m/{m)'^. Indeed, Eq. ( p4|) is con- 
sistent with a typical size of order unity, in contrast with 
the growing average size {m) ^ \/R, a reflection of the 
anomalous algebraic behavior of the size distribution be- 
low the cutoff size. 



VI. THE SIZE- VELOCITY DISTRIBUTION 

So far, we have addressed velocity and size distribu- 
tions separately. However, size and velocity are coupled 
in a nontrivial manner, and for example, slower clus- 
ters should be larger than faster ones. We thus consider 
Pm(v,t), the distribution of clusters of size m and veloc- 
ity V. This joint distribution evolves according to 



dPrnjv.t) 

dt 



= R-^[77lPm+l{v,t) - {m-l)Pra{vA)] 

+ R-^S„,,i[Po{v) - P{v, t)] - c{t)P^{v,t) 
+ I dv' J2 P^{v\t)PJ{v,t)■ (42) 

i+j—m 



The car and cluster velocity distributions are sim- 
ply the zeroth and first moment of the size distribu- 
tion, P{v,t) = Mo{v,t) and G{v,t) = Mi{v,t), with 
Ma{v,t) = ^jnrn'^Pm{v,t). Consequently, the respec- 
tive evolution equations are recovered by summation of 
Eqs. (p^ ) over m. Furthermore, integration over the ve- 
locities gives the size distribution and Eq. (|2^) is recov- 
ered by using Prn{t) = J2m Pm{v, t) . 

It proves useful to introduce the auxiliary functions 
Qm{v,t) = dv'Pm{v' ,t). The cluster size distri- 
bution can be expressed through these auxiliary func- 
tions, Pm{t) — Qm{0,t). Additionally, the identity 
Q{v,t)+ J2Qm,{v,t) — R^^ + c{t) holds. Integrating 
Eqs. ( ^4 ) over v gives 

= R-'[mQm+iiv, t) - (m - 1)Q„(«, t)] 

ot 



+ R-^5,n,iq{v, t) ~ c{t)Q,n{v, t) (43) 
"i{v,t)Qj{v,t) 



-2 ^ ' 

i+j—m 



with q{v,t) — dv'[Po{v' ,t) — P{v',t)] or alternatively 
q{v,t) = 1 - J{t) + Q{v,t) - Qo{v). In deriving (gl) 
we used the following boundary conditions: Q„i = 0, 
Qo = 1 + R~^, and Q = c{t) + R~^ at v = 00. Since 
the velocity plays the role of a parameter, Eqs. (^) can 
be treated as ordinary differential equations. We again 
restrict our attention to the steady state where 

CooQni{v) = R^^[mQ„i+l{v) - (to - l)Qm{v)] 

+ R-'s„,My) + \ E Q-(^)QA^)^ (44) 

i+j—m 

with q{v) = qao{v) ^ 1 - Coo + Qoo(w) - Qo{v)- Introduc- 
ing the generating function 



Q(z,i;) = c-i Ez™Q„(«) 



(45) 



reduces Eqs. ( ^ ) into a set (parameterized by v) of Ric- 
cati equations for Q = Q{z, v): 
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2Q + Z 



g(0) + l-coo^ hz\z 



(46) 



This Riccati equation reduces to Eq. (|T]) when v — 0. 
The above treatment of the size distribution suggests 
that the derivative term in Eq. (^) is negligible for suf- 
ficiently small sizes. In this case, Eq. (Eq) simplifies 
to — 2Q + zq{v)/q{0) = which is solved to give 
Q{z,v) = 1 — \/l — zq{v)/q{0). Using the large R be- 
havior q{v)/q{0) ^ 1 + Q{v) — Qo{v) yields 



Q(z, v)^l- ^l-z[l + Qiv)-Qo{v)]. 



(47) 



Expanding the expression on the right-hand side in pow- 
ers of z we arrive at 



,{v) ^ P^[l + Qiv) ~ Qo{v)Y 



(48) 



with P,n the size distribution (|32|). The cluster size- 
velocity distribution is obtained by differentiating the 
auxiliary distribution Qm{v) 



P^{v) ~ mP^[Po{v)-P{v)][l + Q{v)~Qoiv)Y 



(49) 



Similar to the velocity distribution and the relaxation 
scales, the size velocity distribution as well can be ob- 
tained explicitly from the auxiliary function Q{v). Con- 
sequently, it is characterized by a boundary layer struc- 
ture. The size-velocity distribution is characterized by 
an exponential dependence upon the size, with a velocity 
dependent prefactor. Additionally, there is an algebraic 
prefactor that characterizes the overall size distribution. 

The detailed analysis of the cluster size distribution 
suggests that these results apply only for sufficiently 
small sizes. Eqs. (|47|)-(49) should hold as long as the 
(dropped) term i?^^mQ,„+i(w, t) — (m — l)Q„i{v,t)] is 
negligible compared with the (kept) term CooQm- Using 
(Eq) , the above approximation is valid when 



(50) 



m< Vi? [Qo{v) - Q{v)] . 



Hence, the range of validity of Eq. (^) strongly de- 
pends on the cluster velocity. This can be seen us- 
ing the average cluster size {m{v)) ~ G{v)/P{v) = 
Y.m'nT'Pm{v)/Y.mPm{v), given by Eq. (|2|). Estimat- 
ing the same quantity from Eq. (]49|), gives the correct 
leading large R behavior when f 3> u*, while it gives 
a diverging average size rather than {m{v)) ^ R when 
V ^ 0. Indeed, the condition ( [50| ) is satisfied by the 
{m{v)) only outside the boundary layer. Therefore the 
approximate cluster size- velocity distribution is useful for 
small and average sizes when v ^ v*, while it holds only 
for sufficiently small sizes when u <C w*. Obtaining the 
large size tail requires a more detailed analysis similar to 
that performed for the size distribution. 



VII. SUMMARY AND OUTLOOK 

In summary, we introduced an approximation method 
for analyzing the Boltzmann equations for one dimen- 
sional traffic flows. In analogy with the Maxwell model 
(MM) of kinetic theory, we assumed a constant collision 
rate. This approach results in first order (in the veloc- 
ity) differential equations. Analysis of these equations 
leads to explicit expressions for time dependent velocity 
distributions. Size- velocity distributions can be deter- 
mined in the steady state as well. Although there are 
some quantitative deviations, the overall qualitative be- 
havior including a boundary layer structure, existence of 
laminar and congested phases, etc. is in agreement with 
the results of the original Boltzmann equations. Several 
quantities such as the size growth exponent 1/2 actu- 
ally agrees with the Boltzmann equation. We conclude 
that overall, the Maxwell approach is faithful to nature 
of the problem and thus provides a useful approximation 
scheme. 

The MM allows explicit calculation of several impor- 
tant features, which are otherwise difficult to obtain. The 
approach towards the steady state is generally exponen- 
tial and is characterized by a wide spectrum of velocity 
dependent relaxation scales, the smallest of which corre- 
sponds to the overall cluster density. The steady state 
size distribution exhibits an unusual scaling form with 
a scaling variable m/(m)^ rather than m/{m) which is 
naively expected. Additionally, the typical size which 
is of order unity is much smaller than the average size 
which grows with the collision number. This is a conse- 
quence of the algebraically diverging distribution of small 
sizes. This is an outcome of the non-equilibrium nature of 
the steady state that does not satisfy detailed balance as 
passing events reduce the cluster size by one only, while 
clustering events can increase the cluster size by a large 
number. This feature is independent of the details of 
the collision mechanism and we expect the most features 
underlying the size distribution to hold generally. 

The MM can be refined and systematically improved. 
Some of the quantitative disagreements between the 
Maxwell and Boltzmann equation are rather obvious. For 
example the correct value of the crossover velocity can 
be obtained by replacing the integral dv'Po{v') with 
the integral dv'{v — v')Po{v'). This compensates for 
the approximate kernel taken in the MM and results in 
the correct scaling exponents for the crossover velocity in 
both passing and no passing zones. 

Furthermore, an appropriate choice of the value of 
the prefactor uq reduces the discrepancies between the 
two approaches. For example, the MM gives a univer- 
sal dependence of the density upon the collision number, 
c ~ R~-^/^. However, for the BE if one assumes an al- 
gebraic intrinsic distribution near the origin, Po{v) ~ 
as V —^ different behaviors are found for positive and 
negative /x [Ea. For /i > 0, the density exhibits the uni- 



versal behavior, Co 



R, 



-1/2 







, while for /i < the den- 
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sity becomes /i dependent, Co 



-(Ai+l)/(M+2) 



Rq = CoVqIq is the collision number within the Boltz- 



mann framework. Choosing uq 



R 



m/(m+2) 
■0 



actual BE behavior) implies R — vRq ^ _Rq 



(the 

(2m+2)/(m+2) 



and thence Co 



R: 



-(M+l)/(Ai+2) 



R-^/"^. Therefore the 



BE and MM results are consistent with each other if the 
appropriate choice for the collision rate uq — (v) is made. 

Additionally, it would be interesting to compare the 
MM with the actual traffic process. Although the BE 
description is plausible at the steady state, it is clearly 
an approximation for the transient regime. For exam- 
ple, the BE differs from the exact behavior in the no 
passing case. Another avenue for further research is 
inhomogeneous traffic flows where a hydrodynamic de- 
scription may prove useful. The hydrodynamic frame- 
work should involve a multicomponent fluid parametrized 
by the cluster size m. Specifically, the macroscopic 
description requires the density Pm(x,t), the average 
velocity Um(x,t) = P^~^{x,t) J dvvPm(v,x,t), and the 
"temperature" (the average velocity square) for each m. 
Such infinite-fluid hydrodynamics may be quite different 
from the conventional one-fluid hydrodynamics. Indeed, 
Eq. ( ^sl) shows that the velocity decreases as the mass 
increases. Similar results apply for the temperature and 
thus, the equipartition of "energy" breaks down as well, 
in contrast with usual hydrodynamics. 

In conclusion, the MM is a useful approximation to the 
kinetic traffic equations. This approach may be applica- 
ble to other traffic problems as well. In particular, it will 
be interesting to apply the Maxwell approach to traffic 
models with more realistic passing mechanisms. 

The work of PLK has been supported by NSF (grant 
DMR-9632059) and ARO (grant DAAH04-96-1-0114). 
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APPENDIX A: THE CONDITIONAL VELOCITY DISTRIBUTION 



The car velocity distribution involves the leading car as well as the rest of slowed down cars in the cluster. The 
former is described by the cluster velocity distribution, while the latter is represented by P{v,v',t), the density of 
cars of intrinsic velocity v driving with velocity v' . For the Maxwell model, the master equation for this conditional 
distribution reads 



R-^P{v,v',t) + Piv,t)P{v',t) + P{v',t) dv"P{v,v" ,t) - P{v,v' ,t) I dv"P{v\t). (Al) 



dP{v,v',t) 
di 

The first term accounts for loss due to escape, while the rest of the terms represent changes due to collisions. Integrating 
these equations over the first velocity index and using the relation G(v) = P{v) + dw P{w, v), one indeed recovers 
the rate equation (p^). 

Let us introduce the auxiliary function Q{v, v' , t) — J^, dw P{v, w) which gives the conditional velocity distribution 
by differentiation P{v,v',t) = —dQ{v,v',t)/dv'. This auxiliary function evolves according to 



Integrating Eq. ( A2) over v' and using the boundary condition Q{v, v, t) = 0, we get 

dQ[v,v',t) 



dt 



= Q{v', t)Q{v, v', t) + P{v, t)Q{v', t) - P{v, t)Q{v, t). 



(A2) 



(A3) 



This is a linear inhomogeneous differential equ ation for the auxiliary function Q(v,v' ,t) which includes already known 
cluster velocity distributions. Integrating Eq. (|A3|) we arrive at 



Q{v,v',t) 



dt' P{v,t')[Q{v,t')~Q{v',t')] exp 



dt"Q{v',t") 



(A4) 



The exact solution (A4) can in principle be reduced to a mo re ex plicit expr essi on by following the procedure detailed 
in Appendix B for transforming the formal solution of Eq. (Bl) into Eq. (B2). Such a solution is very cumbersome 

so we do not give it here. 

The steady state conditional distribution is obtained immediately from Eq. ( [A^ ) Q{v,v') — P{v)[Q{v)/Q{v') — 1]. 
The joint distribution is found by differentiation, P{v,v') — —dQ{v,v')/dv' — P{v)P{v')Q{v) /Q'^{v'), or explicitly 



P{v,v') = 



RPo{v)Po{v' ) 
[l + 2RIo{v')] 



3/2 ■ 



(A5) 



In the laminar phase, this conditional distribution is proportional to R, while it is algebraically suppressed in the 
congested phase. One can verify that Eq. ( [A5| ) is consistent with the cluster and car distributions using the relations 
P{v) — Poiv) — /g dv'P{v,v'), and G{v) — P{v) + J^dw P{w,v), respectively. 



APPENDIX B: THE AUXILIARY CAR VELOCITY DISTRIBUTION 



The master equation (M) for g{v,t) can be integrated formally 



g{v,t) = go{v) 



exp^-^ dt' Q{v,t')j + R-^ dt' exp dt" Q{v, 



t") 



(Bl) 



To obtain more explicit results, we notice that the velocity plays the role of a parameter in Eq. (Bl). We thus change 
variable from t' to q= Q{v, t') and for example 



dt'Q{v,t') 



dq ■ 



Qa ' 1^ ~ Qlo 



In 



This change of variables allows us to perform the integration 
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9{v,t) = go 



= go 



= 90 



Qo 



Ql 



Q'-Q 



R 



Q 



dq 



n- 

Q 



+ 



Q 



,0 

Qo Q^- 



n2 



+ 



RQlo 



RQl 
Q + {Q^-Ql 



RQlo Ql 
Qo 



In 



1 

To 



RQ 
t 
2 



+ 



Q^-Q 
2RQ 



2RQ 
' ln(e-*« 







Qcc iQo ^ Qc 



+ Q c 



(B2) 



In the above derivation we used the identities goo = 9o/RQooj Qo ~ Qoo — ^O' ^^'-^ ~ Qo/RQ^] = lo/RQ 



2 

oo* 
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